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ABSTRACT 


We propose an in-situ formation model for inverse-polarity solar prominences 
and demonstrate it nsing self-consistent 2.5-dimensional magnetohydrodynamics 
simnlations, including thermal conduction along magnetic helds and optically 
thin radiative cooling. The model enables us to form cool dense plasma clouds 
inside a flux rope by radiative condensation, which is regarded as an inverse- 
polarity prominence. Radiative condensation is triggered by changes in the mag¬ 
netic topology, i.e., formation of the flux rope from the sheared arcade field, 
and by thermal imbalance due to the dense plasma trapped inside the flux rope. 
The flux rope is created by imposing converging and shearing motion on the ar¬ 
cade held. Either when the footpoint motion is in the anti-shearing direction or 
when heating is proportional to local density, the thermal state inside the hux 
rope becomes cooling-dominant, leading to radiative condensation. By control¬ 
ling the temperature of condensation, we investigate the relationship between 
the temperature and density of prominences and derive a scaling formula for 
this relationship. This formula suggests that the proposed model reproduces the 
observed density of prominences, which is 10-100 times larger than the coronal 
density. Moreover, the time evolution of the extreme ultraviolet emission syn¬ 
thesized by combining our simulation results with the response function of the 
Solar Dynamics Observatory Atmospheric Imaging Assembly hlters agrees with 
the observed temporal and spatial intensity shift among multi-wavelength during 
in-situ condensation. 


Subject headings: Sun: filaments, prominences 
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1. Introduction 

Solar prominences are cool dense plasma clouds that form in the hot tenuous corona. 
Dense plasmas are sustained at the dipped (concave-up) structures of coronal arcade fields 


(Kippenhahn & Schliiter 1957; Kuperus & Raadu 1974). Prominences can be classified 


into two groups — normal- and inverse-polarity prominences — depending on their 
interior magnetic polarity with respect to the overlying arcade field. In normal-polarity 
prominences, the magnetic polarity matches that of the overlying arcade held. Normal- 
polarity prominences are typically modeled by the Kippenhahn-Schliiter (KS) model. In 
inverse-polarity prominences, the interior magnetic polarity opposes that of the overlying 
arcade held. This prominence group is described by the Kuperus-Raadu (KR) model, 
in which the lower half of the hux rope indicates the inverse-polarity. Both normal- and 


inverse-polarity prominences are observed in the solar atmosphere (Leroy et ah 


An unsolved problem regarding prominences is the origin of their cool dense plasmas. 
Though several models have been proposed, the actual mechanism by which prominences 


form remains debatable (Mackay et ah 2010). One candidate mechanism for prominence 


formation is radiative condensation (Antiochos et ah 1999). Several models for triggering 


radiative condensation have been proposed and demonstrated by numerical simulations. 


One is the evaporation-condensation process (Mackay et ah 2010), in which evaporated 


flows from the chromosphere lead to radiative condensation in the coronal loops. The 
mechanism and constraints of this model have been investigated using one-dimensional 


hydrodynamics simulations (Karpen et ah 2006 Luna et ah 2012) and two-dimensional 


self-consistent magnetohydrodynamics (MHD) simulations (Xia et ah 2012). These studies 


have shown the formation of normal-polarity prominences categorized by the KS model. 


whereas Keppens & Xia (2014) showed that the evaporation-condensation model eventually 


evolved into the funnel prominence in which the flux rope structure is embedded. 
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Choe & Lee (1992) demonstrated the formation of a normal-polarity prominence by a 


different process. In their two-dimensional MHD simulations, the shearing motion at the 
footpoints of the potential arcade fields leads to radiative condensation due to the expansion 


of the arcade fields. Linker et ah (2001) demonstrated the formation of inverse-polarity 


prominences. In their simulations, the chromospheric plasmas are levitated as a result of 
flux rope formation by the magnetic cancellation and condensation occurs. 

Condensation models triggered by evaporation or levitation rely on the injection 
of chromospheric plasmas. On the other hand, a recent observation of a prominence 
formation event showed no clear evidence of chromospheric plasma injection, and claimed 


the presence of in-situ condensation (Berger et al. 2012). Xia et ah (2014a) modeled in-situ 


condensation inside a flux rope by three-dimensional simulation, and showed that the 
features of prominences are well reproduced. Because they used the mechanically stable flux 


rope model obtained from the isothermal simulation (Xia et al. 2014b) and started their 


simulation with an ad hoc thermal imbalance, the process required to trigger the radiative 
condensation is still unclear. 


In this study, we model the process from flux rope formation to radiative condensation 
using self-consistent 2.5-dimensional MHD simulations, including thermal conduction and 
radiative cooling. By comparing the results of multiple simulations with different settings 
regarding the formation model of flux rope and the coronal heating model, we discuss the 
condition for triggering radiative condensation and the possible relationship between the 
temperature and density of the prominence condensations. 

Section 2 introduces our prominence formation model, and Section 3 describes the 
numerical simulation setup. The simulation results are presented in Section 4. The 
condition for triggering radiative condensation and a scaling formula from the possible 
temperature-density relationships are discussed in Section 5. 
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2. Possible processes for triggering in-situ condensation 


The schematic picture of our model for triggering in-situ radiative condensation is 
presented as follows. Initially, the coronal arcade exists in thermal equilibrium (Fig. UKa)). 
The arcade held is transformed into a hux rope structure by converging and shearing motion 


(van Ballegooijen & Martens 1989; Martens & Zwaan 2001). The relatively dense plasmas 


in the lower corona due to stratihcation are trapped inside closed loops of the hux rope, 
where they are elevated into the upper corona (Fig. (b)). The dense plasmas increase the 
radiative cooling in the hux rope, leading to thermal imbalance. Because the closed held 
lines are disconnected from the exterior region, the heat loss by radiative cooling inside the 
hux rope is not conducted outward, and radiative condensation can be triggered. (Fig. 

(c)). 


coronal arcade 


flux rope 
formation 


low density 
stratified 
high density 


1 



(a) Initial state 




(c) Condensation 


Fig. 1.— Schematic of a possible process of in-situ radiative condensation 


3. Numerical Settings 

To demonstrate the proposed model, we perform 2.5-dimensional MHD simulations. 


including thermal conduction along the magnetic held lines, radiative cooling, and gravity 
in a Cartesian {x,y) domain. To investigate the necessary condition for triggering radiative 
condensation, we test diherent direction of the converging and shearing motion for the 
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formation of a flux rope and different coronal heating models. 

3.1. Initial and boundary conditions 

The corona model is set up in a rectangular box, in which the x- and y-axes 
are horizontal and vertical, respectively, and the ; 2 -axis is orthogonal to the x-y 
plane. The corona is stratihed under uniform temperature (T^or = IMK) and gravity 
(fi'cor = 2.7 X 10"^ cm s“^) conditions, 

p{y) = Pcorexp - y , (1) 

Kb J- cor 

p{y) = —p{y)Tcor, (2) 

m 

where pcor = 3.2 x 10“^® g cm“^, ks is Boltzmann’s constant, and m is the mean molecular 
mass. We set m = nip with nip being the proton mass. The force-free arcade field iss 
described as 



where 5^ = 3 G is the held strength at the footpoint. La = 12 Mm is the width, and 
a = 31 Mm is the magnetic scale height of the arcade held. Initially, the system exists in 
mechanical equilibrium. 

The left and right boundaries are subjected to symmetric (for p,p,Vy,By) or anti¬ 
symmetric (for Vx,Vz, Bx, Bz) boundary conditions. A free boundary condition is applied to 
the top. In the region below y = 0, the converging and shearing motions are introduced. We 
test three types of the footpoint motions. One is the converging motion without shearing; 
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the others are converging motion with shearing that increases the magnetic shear of the 
arcade held or with anti-shearing that decreases the shear of the arcade held. In all the 
cases, the velocity components Vx and Vy within this region are set as follows, 


Vx = -vo{t)-^, Vy = 0, {0 < X < La/4:, y <0) (6) 

Vx = Vy = 0, {La/4:<x <La, y <0), (7) 

' noo, (0 < t < ti) (8) 

Voit) = < Vooit2 - t)/{t2 - ti), {ti<t < t2) (9) 

.0, {t>t2) (10) 


where noo = 12 km/s, ti = 1000 s, and t 2 = 1500 s. In the case of no shearing motion, 

Vz = 0. The shearing and the anti-shearing motions are set as follows, 

Vz = ±vo{t)-^, (0 < a; < La/4, 1/< 0) (11) 

-^a/4 

Vz = ±t^o(t) ^g2 (La/4 < x < La, 1/< 0) (12) 

where the plus and minus signs represent the shearing and anti-shearing motions, 
respectively. The magnetic helds are computed with the induction equation by coupling 
with the given converging and shearing motions. The magnetic helds at the bottom 
boundary are hxed to the initial state. The gas pressure and density are hxed by assuming 
hydrostatic equilibrium at a constant temperature of T^or = 10® K. 


3.2. Basic equations 


The basic equations are 


•^eth 


-F V ■ Veth 


— + vVp= -pV ■ u, 

-(eth + p) V ■ u + V ■ {nT^/^hh ■ VT) - L, 


(13) 


dt 


(14) 



Cth — 


p 


7-1’ 
1 , 


dv 1 

— + v ■ Vv = —Vp ^ 
ot p Anp 


V X B) X B + g, 


dB „ ^ 

-^-cVxE, 


E = —V X B 

c 


Airp 


J = —V X B, 
Att 


(15) 

(16) 

( 17 ) 

(18) 
(19) 


where k = 2 x 10“® erg cm“^ s“^ is the coefficient of thermal conduction, h is a unit 

vector along the magnetic held, and y = (0, —pcor, 0) is the gravitational acceleration and p 
is the magnetic diffusion rate. The temperature is computed by the following equation of 
state, 

m r) 

( 20 ) 


^ m p 

T = - -. 

ksp 


For fast reconnection, we adopt the following form of the anomalous resistivity (e.g. 


Yokoyama & Shibata 1994), 


0 , 


(J < Je) 


p = 


( 21 ) 

70 {J/Jc - 1)' , {J > Jc) (22) 

where po = 3.6 x 10^^ cm^ s“^ and Jc = 25 dyn^/^ cm s“^. We restrict p to Pmax = 


18.0 X lO^'^ cm^ s L in Eq. (14) represents the net cooling term, which is expressed as 


L = n^A{T)- H -pJ^, 


(23) 


where n = p/m is the number density, A(T) is the radiative loss function of optically thin 
plasma, H is the background heating rate, and pJ^ represents ohmic heating. We adopt a 
simplihed radiative loss function in Fig. (a) (Hildner 1974) and simulate two different 
coronal heating models. In one model, the heating rate H depends on the local magnetic 
energy density (magnetic pressure) at temperatures above the cutoff temperature Tc and is 












balanced out by the cooling rate when it falls below 

a^Pm, (T > Tc) 
n^AiT), (T<Te) 
Initial thermal equilibrium required that 


H = 


A(Tcor) 

/Q - 

Bl Sir 


mgc 


^bTcov 


(24) 

(25) 

(26) 


where Ucor = Pcor/'^ = 2.0 x 10® cm Selecting a = ksTcor/irngcov), we obtain 


H = 


am = 6.2 X lO”'^ cm^ (i.e., a constant). Equation (25) prevents the temperature from 
decreasing below Tc. In the other model, the heating rate H depends on the local density 
and height at temperatures above Tc and balances with the cooling rate below the cutoff 
temperature Tc. 

ad{y)n, (T > Tc) (27) 

n‘^A{T). (T<Tc) (28) 

For the initial thermal equilibrium, we have ad{y) = ncorA{Tcor) exp[—mgcory/(kBTcor)] 
where ?7,corA(Tcor) = 1.1 x 10“^^ erg s“^. Note that ohmic heating is always much smaller 
than background heating in these settings. Tc is one of the model parameters. For each 
heating and shearing model, we vary Tc as 0.3 MK, 0.4 MK, and 0.5 MK. The investigated 
cases are summarized in Table [U 


3.3. Numerical scheme and grid size 


The numerical scheme is based on the CIP-MOCCT method (Kudoh et ah 1999). To 


prevent numerical oscillations in the advection term (left-hand side of Eqs. (13)-(16)), we 


apply the rational CIP scheme (Xiao et ah 1996), in which the rational function replaces 


the cubic function as an interpolation function of CIP. The thermal conduction, radiative 
cooling, and heating terms are included as source terms in non-advection phases of CIP 
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method. Thermal conduction is explicitly solved by a second-order accurate slope-limiting 
method (Sharma & Hammett [2007 ). 


The required grid size for preventing numerical thermal instability can be estimated by 


the formula presented in Koyama & Inutsuka (2004). When applying this formula to our 
simulation, we assume that the density of the condensations is 10 times the coronal density, 


A,. _ 1 / kTV'^ 

3 “3Vn2A(T,)’ 


(29) 


where A is the grid size, Xp is the Field length, Tc is the cutoff temperature (also assume 
to be the temperature of condensations), and ric = lOucor is the assumed density of 
condensations. Figure]^ (b) shows the required grid size against temperature. Under the 
constraint of Eq. ( [2^ , the grid size is determined to be A = 120 km, 60 km, and 30 km 
in the cases of Tc =0.5 MK, 0.4 MK, and 0.3 MK, respectively. The horizontal grid size 
Ax is uniform. The vertical grid size Ay is uniform up to ?/ < 24 Mm and increases to 
Ay = 1.8 X 10^ km ai y > 24 Mm to ensure a sufficiently distant upper boundary. 
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Fig. 2.— Simplified radiative cooling functions and required grid sizes 
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Results 


The right-most column of Table [T] shows our results in the presence of the radiative 
condensation. We hnd that the necessary condition for the radiative condensation is either 
anti-shearing motion or the heating proportional to the local density. Panels (a)~(d) of 
Fig. [^are snapshots of the time evolution and Fig. |^is the hnal state in case Ml, which 
is a the typical case for the radiative condensation. The initial state exists in mechanical 
and thermal equilibrium (Fig. [^(a)). Converging motion triggers reconnection above the 
PIL at a; = 0, and a flux rope is formed (Fig. (b)). As the reconnection proceeds, the 
relatively dense plasmas in the lower corona are trapped inside the flux rope and lifted to 
the upper corona. The dense plasmas enhances radiative cooling in the flux rope, which 
overwhelms the background heating and leads to thermal imbalance (Fig. (c)). Because 
thermal conduction works only along closed magnetic held lines of the hux rope, it can not 
compensate for the thermal imbalance. Consequently, radiative condensation is triggered, 
and the cool dense plasmas accumulates in the lower part of the hux rope (Fig. 1(d)). 
Figure (b) shows prohles of temperature, number density, and horizontal component of 
the magnetic held along the y-axis in the quasi-static state (Fig. (a)). The temperature 
of the prominence corresponding to Tc is approximately 0.3 MK in this case. The heating 


settings (Eqs. (25) and (28)) maintains the temperature of the cool dense plasmas at around 
Tc. According to Fig. (b), the magnetic polarity of the cool dense region opposes that 
of the initial arcade held. Figures (a) and (b) clearly show that radiative condensation 
eventually establishes inverse-polarity prominence in the hux rope. Also, the high-density 
region (warm color region in the density plot of Fig. (a)) resembles the tower prominence 


(Berger et al. 2012), and the low-density region surrounding the prominence (dark region 


in the same plot) can be regarded as the coronal cavity. Figure shows each force at a 
time of =4250s. The mechanical quasi-static state is realized in the vertical direction. 
The upward magnetic tension force is balanced with the downward pressure gradient force 
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and gravity, i.e., the prominence mass is snpported mainly by the magnetic tension. This 
force balance system is also fonnd in the parametric stndy by adiabatic MHD simulations 


of Hillier & van Ballegooijen (2013). Our result mostly corresponds to their case (e) of 
model 1 regarding the force balance, whereas the prominence density and plasma (3 in our 


simulation are lower and higher than those in Hillier & van Ballegooijen (2013), respectively. 


They revealed that the magnetic tension force is given by the vertical stretching of the flux 
rope, and the curvature of the held lines inside the prominence can be smaller because 
the sufficiently large magnetic tension comes from the larger held strength by compression 
of the prominence mass. Our simulation realizes these features of prominence support as 
shown in Fig. |^(a). Figures]^ shows the time evolution of case M5, which is a typical case 
with no radiative condensation. The hux rope is heated-up rather than cooled-down in this 


case. 


5. Discussion 

As shown in Table the radiative condensation occurs either when the anti-shearing 
motion is imposed or when the heating model is proportional to the density, i.e., H = adp. 
The radiative condensation is triggered by the cooling-dominant thermal imbalance in the 
flux rope composed of the thermally isolated closed loops. The interior of the flux rope 
is relatively dense because of the converging motion as well as the initial stratihcation. 
Because the radiative cooling inside the flux rope is enhanced by the dense plasma, the 
cooling-dominant thermal imbalance is reasonable if the heating changes only slightly. 
Because we set the heating depending on the magnetic pressure (cases M1-M5, hereafter 
called M-cases) and local density (cases D1-D5, hereafter called D-cases), the change in 
heating rate is not negligible. In some cases (M5 and M6), a heating-dominant thermal 
imbalance is created and no condensations occur. In the following subsections, we discuss 
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case 

Heating 

Shearing 

Te (MK) 

Condensation 

Ml 


— 

0.3 

Yes 

M2 


— 

0.4 

Yes 

M3 

H cXfYiPffi 

— 

0.5 

Yes 

M4 


0 

0.5 

No 

M5 


+ 

0.5 

No 

D1 


— 

0.3 

Yes 

D2 


— 

0.4 

Yes 

D3 

H = adP 

— 

0.5 

Yes 

D4 


0 

0.5 

Yes 

D5 


+ 

0.5 

Yes 


Table 1: The presence of radiative condensation in each case. The second column shows the 
heating model. The third column shows the shearing model; plus sign (+), minus sign (—), 
and 0 represent the shearing, anti-shearing, and no shearing cases, respectively. The fourth 
column shows the cutoff temperature. The fifth column shows the presence of condensation. 
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Fig. 3.— Panels (a)-(d) are snapshots of the time evolution in case Ml. Colors represent 
temperature and density, white lines represent magnetic held lines, and arrows represent the 
velocity helds. 
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(a) 


time= 4250s 




6 

X (Mm) 


6 

X (Mm) 


0.2 0.8 1.4 8.8 9.3 


9.8 

.■3i 



temperature (MK) log.|Q[n cm ’ 

Fig. 4.— Panel (a) is the final state in case Ml. Colors, white lines and arrows represent the 
same quantities as in Figj^ Panel (b) shows the profiles of temperature, number density, 
and Bx along the ?/-axis. The solid and dashed lines are the profiles in the final and initial 
states, respectively. 


the mechanism for creating the cooling/heating-dominant thermal imbalance in the M- and 
D-cases. Then, our model is compared with previous theoretical studies and observations. 


5.1. Thermal imbalance in M-cases 

In M-cases where the heating rate is proportional to the magnetic pressure, the 
anti-shearing and shearing motion causes a larger discrepancy between the cooling 
and heating rates, establishing the cooling- and heating-dominant thermal imbalances, 
respectively. Figure shows the cooling and heating rates along the |/-axis at time=500s in 
case M3-M5. In case M3 (anti-shearing case), the cooling rate is the highest of the three 
cases, whereas the heating rate is the lowest. Case M5 (shearing case) has the opposite 
tendency, where the cooling rate is the lowest and the heating rate is the highest. We 
discuss how these tendencies are established from a perspective of the evolution of the 
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0 5 10 15 20 25 

y (Mm) 


Fig. 5.— The vertical components of the forces along the y-axis. The dotted line, dashed line, 
and dot-dashed line represent the pressure gradient (the sum of the gas and magnetic pressure 
gradients), magnetic tension force, and gravity, respectively. The solid line represents sum 
of all the forces. 
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25 
20 
15 
10 

5 
0 

0 6 12 
X (Mm) 


time= 2500s 



0.2 0.8 1.4 8.7 9.2 9.7 0.2 0.8 1.4 8.7 9.2 9.7 

temperature (MK) log.|Q[n cm'^] temperature (MK) log.|g[n cm'^] 

Fig. 6.— Panels (a)-(b) are snapshots of the time evolntion in case M5. Colors and arrows 
denote the same qnantities as those in Figj^ Black lines represent magnetic held lines. 


magnetic held and its inhnence on the cooling and heating rates. 


In the anti-shearing case, the hnx rope is pnshed npward by the reconnection and 
conhned by the downward magnetic tension of the overlying arcade held. Fignre shows 
the vertical tension force along the y-axis at time=500 s in case M3 (anti-shearing), M4 (no 
shearing), and M5 (shearing). The downward magnetic tension force in the anti-shearing 
case is larger than those in the shearing and no shearing cases. The reason the anti-shearing 
case has the larger downward tension force is as follows: the shear angle is related to the 
height of the arcade held as 


a = 


71 COS 6 ’ 

where a is the arcade height, 6 = aTcta.n{Bz/B^) is the shear angle against the positive 
x-axis in the x^-plane. La is the arcade width, and we assnme a linear force-free arcade 


(30) 


held for simplicity. Eqnation (30) indicates that the redaction of magnetic shear by the 


anti-shearing motion makes the arcade held shorter. Dne to the shrink of the arcade held, 
the hnx rope in the anti-shearing case experiences a larger downward magnetic tension 
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force. 

The flux rope in the anti-shearing case is pinched by the overlying arcade held and 
pushed by the upward magnetic tension of the reconnected loops, leading to a density 
enhancement by compression. Figure shows the time evolution of area and average 
density inside the dash-dotted line in Fig. The average density is computed by dividing 
the total mass inside the dash-dotted line by the corresponding area. The total mass is 
almost the same between these three cases and is conserved after the hux rope is detached 
from the bottom boundary. In Fig. (b), the anti-shearing case shows the enhancement of 
the average density. This result indicates that the density enhancement inside the flux rope 
in the anti-shearing case is due to the compression by the magnetic tension as well as the 
converging motion. 

The anti-shearing motion lead to the higher cooling rate as well as the lower heating 
rate. The anti-shearing motion reduces the magnetic component parallel to the PIL, 
corresponding to in our simulations, and decreases the magnetic pressure, leading to the 
lowest heating rate. Thus, the interior of the flux rope in case M3 suffers from unstoppable 
cooling most efficiently. In case M5 (shearing case), on the other hand, the heating rate 
increases as the magnetic pressure increases, and the enhancement of the cooling rate 
is too small to overwhelm the heating rate. Eventually, the shearing motion is likely to 
establish the heating-dominant thermal imbalance, resulting only in the formation of hot 
flux rope. These results suggest that when the coronal heating depends on the magnetic 
energy density, the continuous shearing motion can not trigger the radiative condensation, 
whereas the temporary anti-shearing motion can indeed trigger it, leading to the prominence 
formation. 
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5.2. Thermal imbalance in D-cases 

From the results of the D-cases in which the heating is proportional to the density, it is 
found that the cooling-dominant thermal imbalance can be achieved only by the difference 
between the heating and cooling rates in the dependence on the density. Figure [T^ shows the 
cooling and the heating rates in the cases D3-D5. The cooling rate in case D3 (anti-shearing 
case) is highest because the density enhancement is largest in the anti-shearing case, as 
discussed in the previous subsection. Unlike case M3, however, the heating rate in case 
D3 is also highest. The combination of the higher cooling rate and the higher heating 
rate can be a disadvantage to enhancing the thermal imbalance. What actually sets the 
cooling-dominant thermal imbalance in D-cases is the different density dependence between 
the heating {H oc p) and cooling terms {R oc p^). The increment of the heating rate against 
the density increase is always smaller than that of the cooling rate. Consequently, the 
cooling overwhelms the heating, and the radiative condensation is triggered in all D-cases. 
These results suggest that if the coronal heating is only proportional to the density, the flux 
ropes in the corona are likely to cool down. The exact modeling of the coronal heating is 
required to conhrm this mechanism. 


5.3. Comparison with other theoretical studies 


In evaporation-condensation models (Antiochos et ah 1999 Karpen et al. 2006), 


radiative condensation is triggered by the transport of dense plasma from lower altitudes. 
Our model is consistent with this mechanism, but it adopts a different approach to suppress 
the thermal conduction effects. In the evaporation-condensation model, thermal conduction 
effects are suppressed by preparing long magnetic field lines, whereas in our model, they 
are suppressed by closed loop formations. In other words, our model altered the magnetic 


held topology. Both our study and that of Choe & Lee (1992) assume that photospheric 
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motion triggers the prominence formation. However, the direction of motion differs between 


the studies. In Choe & Lee (1992), radiative condensation arises by shearing motion 


along the PIL and subsequent expansion of the arcade field. In our study, converging 
motion towards the PIL and subsequent formation of closed magnetic fields (flux rope) are 
essential precursors of radiative condensation. Moreover, whereas their models reproduce 
normal-polarity prominences, our models form inverse-polarity prominences. The formation 
of the inverse-polarity prominence is found both in the results of our model and that of 


Linker et ah (2001). The difference is that our model does not require chromospheric 


plasma injection to trigger the radiative condensation; hence, our model reproduce the 


in-situ condensation. Xia et ah (2014a) performed three-dimensional simulation of the 


in-situ condensation inside the flux rope. The flux rope system in Xia et ah (2014a) was 


created by the converging and shearing motion in the isothermal simulation of Xia et ah 


(2014b). Because the strategy used to establish cooling-dominant thermal imbalance in Xia 


et ah (2014a) was the parameterized heating, it was still unclear why the thermal imbalance 


was created inside the flux rope. We ensure that the cooling-dominant thermal imbalance 
can be created inside the flux rope by two different mechanisms: the anti-shearing motion 
or the different dependence on the local density between the heating and cooling rates. We 
find that the direction of the shear motion strongly affected the thermal processes when 
the heating depends on the magnetic energy density, and it controls the presence of in-situ 
condensation. When the heating is proportional to the density, on the other hand, the 
cooling-dominant thermal imbalance is always established by the difference between the 
cooling and heating terms in the dependence on density, leading to radiative condensation. 

The formation mechanism through the anti-shearing motion smoothly connects to some 
of the eruptive models. These include instability models such as kink and torus instability 


(Kliem et ah 2004 Torok & Kliem 2005; Kliem & Torok 2006 Fan & Gibson 2007). The 


kink instability is likely to be triggered because the anti-shearing motion increases the twist 
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of the flux rope. The torus instability also £t because the reduction of shear can alleviate 
the critical height of torus instability. The other model that £t with our mechanism is 


the reversed shear model of Kusano et ah (2004). We stopped the anti-shearing motion 
before the shear reversal. If we continued to impose the anti-shearing motion, the flux 
rope formation, radiative condensation, and eruption would be successively occurred even 
without the converging motion. Both the shearing and anti-shearing motions have been 
considered as the possible process of eruptions (Kusano 2002 Amari et al. [2003 ). Our 
results suggest that continuous shearing motion results in the eruption of the hot flux 
rope corresponding to flare eruption, whereas the anti-shearing motion causes the radiative 
condensation of the prominence formation. The eruptive mechanism coupling with flux 


emergence is also the subject to investigate (Chen & Shibata 2000 Shiota et al. 2005 


Kaneko fc Yokoyama||2014). The results of our simulation with the good mechanical balance 


can be directly adopted as the initial condition for these eruptive studies. The model 
including the dense materials of the prominence may reveal the more detailed structure of 
coronal mass ejection such as the dense core and its evolution. 




0 5 10 15 20 25 
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Fig. 7.— Cooling and heating rates along the y-axis at time=500s in Figj^ Solid, dashed, 
and dot-dashed lines represent cases M3, M4 and M5, respectively. 
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Fig. 8.— The vertical magnetic tension force along the y-axis at time=500s. Solid, dashed, 
and dot-dashed lines represent M3, M4, and M5, respectively. 




time (s) time (s) 

Fig. 9.— Time evolution of area and average density inside the dash-dotted line in Figj^ 
Solid, dashed, and dash-dotted lines correspond to M3, M4, and M5, respectively 
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Fig. 10.— Cooling and heating rates along the y-axis at tinie=500s in Figj^ Solid, dashed 
and dot-dashed lines represent cases D3, D4 and D5, respectively. 


5.4. Comparison with observation 


We compare the simulated process of in-situ condensation with the observation by 


Berger et ah (2012) by synthetic emission through the extreme ultraviolet (EUV) hlters of 


the Solar Dynamics Observatory Atmospheric Imaging Assembly (SDO/AIA). The emission 
of a certain wavelength channel i is expressed as 


Di = J nlKi{T)dl, (31) 

where Dj, Ue, Ki{T), and I represent the pixel response to the photon flux, electron number 
density, temperature response function of the AIA filters, and the distance along the 
line-of-site (LOS), respectively. The temperature response functions are obtained from 
aia_get_response.pro in the Solar Software Library. To compute the synthetic emission, we 
choose the z-axis as the LOS. In our 2.5-dimensional simulation, the physical variables are 
assumed to be constant along 2 ;-axis; hence, the integration along the LOS is altered to the 
multiplication of an arbitrary constant Iz as follows. 


A = nlKt(T)k. 


(32) 
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Figure 11 is the time evolution of the synthetic emission of each EUV hlter. Note that Fig. 


11 covers the time only before the temperature reaches the cutoff. Figure 11 clearly shows 


that the areas of strong emission progressively shift from the 193A channel, through 17lA 
channel, to 13lA channel both in time and altitude. The emission shift was claimed to be 


evidence of in-situ condensation in Berger et al. (2012), and our simulation has validated 


their argument. The shift from the 21lA channel to the 193A channel is not so clear 


because the initial coronal temperature is set to IMK. The shift in altitude in Fig. 11 


results from the descent of the flux rope and high density plasma. It can also be realized 


by some other mechanisms such as magnetic reconnection in the weighted dip (Liu et al. 


2012 Keppens & Xia 2014), or cross-held slippage (Low et al. 2012a||b ). These mechanisms 
are not realized in our simulations because we set the cutoff temperature such that the 
simulated prominence density can not be sufficiently large. 

Figure [T^ shows the relationship between the temperature and the maximum density of 
the condensations obtained from our results. The lower the condensation temperature, the 
higher the density of the condensates. The resultant density is independent of the heating 
model. The temperature and density of the prominence are related by the following scaling 
formula. 


n 


= 85 


T 


-1 


. 2 X 109 cm-3 ) V 10^ K 
The scaling relationship (Eq. (p33|) is plotted as a solid line in Fig. 12 This formula 


(33) 


indicates that the dense condensation results from re-strati£cation along the magnetic held 
lines with temperature T^. Because the mass is conserved along the held lines and conhned 
to the length scale of the scale height ksTd{mgcor), the condensate density is inversely 
proportional to temperature. From Eq. ( [M| ), we concluded that our model can reproduce 
the observed prominence density, which is 10-100 times larger than the surrounding coronal 
density at typical temperatures (below 10^ K). 
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The direct simulation without the temperature cutoff is necessary to verify the realistic 
prominence density and compare DEM with observation. Including the chromosphere 
in our model is also future work. Due to the thermal conduction from the corona to 
chromosphere, the background heating rate must be larger than the cooling rate. When the 
initial coronal density is lower, the magnitude of the initial thermal imbalance is lower, and 
the present model without the initial thermal imbalance will be validated. It is the next 
step for our model to be extended to three dimensions. Self-consistent three-dimensional 
MHD simulation including coronal thermal conduction is still a challenging issue. The 


comparison with observations including the LOS effect (Luna et ah 2012; Xia et ah 2014a) 
must contribute to the further understanding of the prominence formation. 


6. Conclusion 

We proposed a new in-situ formation model of inverse-polarity prominence. The model 
was demonstrated in 2.5-dimensional MHD simulations, involving thermal conduction and 
radiative cooling. Either flux rope formation by anti-shearing motion or a heating model 
proportional to the local density is necessary for a in-situ prominence formation. From the 
results of a parameter study on cutoff temperature, we derived an empirical scaling formula 
that related the density of the prominences to their temperature. Our study quantitatively 
demonstrated that our in-situ formation model can reproduce the observed density of 
prominences. The study also explained the observed intensity shift among multi-wavelength 
EUV emission of in-situ condensation. 

We are grateful for the helpful discussion with Dr. C. Xia and Prof. Dr. R. Keppens 
during my short stay in KU Leuven. This trip was supported by the Program for 
Leading Graduate School, MEXT, Japan. This work was supported by JSPS KAKENHI 
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Fig. 11. — Time evolution of synthetic EUV emission through the 21lA, 193A, 17lA, and 
131A channels of AIA. logTpeak = 6.2, 6.1, 5.8, and 5.6, respectively, where Tpeak is the peak 


temperature of each channel. Iz in Eq. (32) is set to 10 Mm. The emission is averaged in 
the area of x < 3Mm. 
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Tc (MK) 


Fig. 12.— Relationship between Tc and maximnm density. The modeled cases are indicated 


by letters. The solid line is the resnlt of applying the scaling formnla Eq. (33). 
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